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^ In [Ij, we have proposed a model where natural inflation is realized on a steep potential 

I (^(0) ~ cos(0//) with / <^ Mp) as a consequence of the interaction of the infiaton with 

2 gauge fields through the couphng F^^, F^^" . In the present work we study the nongaussiani- 

^ ties and the spectrum of tensor modes generated in this scenario. The nongaussianities turn 

, ^, out to be compatible with current observations and can be large enough to be detectable 

by Planck. The non-observation of tensor modes imposes new constraints on the parameter 

^ space of the system that are about one order of magnitude stronger than those found in P] . 

0^ More importantly, in certain regions of the parameter space tensor modes might be detected 

^ by upcoming Cosmic Microwave Background experiments even if inflation occurs at energies 

lO as low as the TeV scale. In this case the tensor modes would be chiral, and would lead to 

cn distinctive parity- violating correlation functions in the CMB. 
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1 Introduction 

An axion-like field $ with 



A' 



4 



COBI- 



(1) 



represents the simplest example of spin-0 degree of freedom with a nontrivial and yet radia- 
tively stable potential. The technically natural flatness of V{^) makes $ an ideal infiaton 
candidate |2]. In all known UV-complete models including potentials of the form (jT|, the 
parameter / is significantly smaller than the reduced Planck mass Mp = I/a/SttG [3]. It is 
actually conjectured [1] that it is impossible for a UV-complete theory that includes gravity 
to fulfill the condition / > Mp: the potential ([I]), while controllably fiat, cannot be arbi- 
trarily fiat. Unfortunately, the condition / > Mp is precisely what ([I]) requires to support 
inflation. 

There exists by now a substantial body of literature (see e.g. [3 El [71 El E]) that aims 
at maintaining the properties of radiative stability of the axion-like field while changing the 
shape of V{^) to allow it to support inflation. In [I], we have taken a different approach. 
Instead of deforming V^($), we have coupled $ to matter. More specifically, we have con- 
sidered the effect on the dynamics of $ induced by its coupling to a f/(l) gauge field A^ 
via 

C^^^=j^F,,F^\ (2) 

with a a dimensionless parameter. In the presence of this coupling, the rolling infiaton uses a 
finite fraction of its kinetic energy to produce quanta of A^. This coupling has thus the effect 
of slowing down the rolling of $, making it possible to achieve inflation even if V{<^) is steep 
(i.e., / ^ Mp). This situation is analogous to that leading to warm inflation [TOl [TTl [T2] . 
to trapped inflation [13] and, more recently, to the mechanism discussed in [13]. We also 
note that a model related the present scenario has been recently proposed in [13 US], where 
the infiaton is coupled to SU{2) gauge fields leading to a classical, rotationally invariant 
configuration able to support slow-roll even for / < Mp. 

The analysis of [1] shows that it is possible to obtain sufficiently long inflation, irrespective 
of the value of /, if a is of the order of a few hundred. In the simplest version of the model 
the amplitude of the scalar perturbations is too large to match observations. This happens 
(almost) irrespective of the parameters in the theory because, once a is large enough to 
allow for sufficient inflation, the observables of this model have a very weak, logarithmic 
dependence on the parameters of the theory (see section 2). It is however possible to reduce 
the amplitude of the scalar perturbations to its observed value by considering a modification 
of the theory where the infiaton couples to A/" ~ 10^ ^(1) gauge fields. 

Is this model in any sense generic? Of course it is not - both the large value of a and the 
large number of gauge fields are not generic predictions of UV-complete theories. However, 
no fully working model of inflation is generic. And actually, in the spirit of chaotic inflation, 
inflation does not need to be generic: as long as some corner of the parameter space of the 
theory allows it, inflation will happen and the inflating part of the multiverse will overwhelm 
the remaining sections. In [T] we have discussed situations where the conditions a = (9(10^) 
and J\f = O{10^) can be consistently realized. In this sense, this set-up is as generic as we 
can require. 



The only way of discriminating between models is to study their predictions and compare 
them to observations. In [T] we have shown that the spectral index of scalar perturbations 
is n ~ 1 + (2/7ra) (/ V" IV'\ and so it is extremely close to 1. More recently, [T71 HE] have 
studied the nongaussianities and [19j has studied the spectrum of (chiral) gravitational waves 
induced by the coupling ([2]), that might even be observable by ground-based gravitational 
interferometers in the near future ^D]. The analyses of refs. [T71 [TSl IIHI EO] (see also [21]) 
are however performed in the regime where the backreaction of the gauge field on the rolling 
inflaton is negligible, and therefore inflation is obtained because V(<l') (that is otherwise 
considered to have an arbitrary shape) is sufficiently fiat. 

In the present paper we study the nongaussianities and the spectrum of tensor modes in 
the model [T]. We find that the j^i^ computed on equilateral configurations takes values in 
the range /^™ ~ —30 -i — 5 depending on the energy scale of infiation. Nongaussianities are 
therefore within the current bounds, and could be detected by Planck. 

Our predictions for the spectrum of tensor modes depend on the value of the coupling a 
and on the scale of infiation. We find that non-observation of tensor modes induces a bound 
a -^ 300 -T- 10^ depending on the scale of infiation. Moreover, as usual, the scale of infiation 
has to be below the 10^^ GeV if we want to avoid overproduction of tensors. Remarkably, 
however, the tensors-to-scalar ratio can be C(10~^), and therefore detectable by upcoming 
CMB polarization experiments, irrespective of the energy scale at which infiation occurs, 
provided a has the correct magnitude. The system of |T] provides therefore an example of 
a model where gravitational waves of infiationary origin could be detected even if infiation 
happens at energies much smaller than the GUT scale [22] • Unlike tensors produced by 
most other mechanisms, these gravitational waves would be fully chiral [IH], a property that 
would show up in parity-violating correlation functions in the CMB, providing a very specific 
signature of this scenario. 

The paper is organized as follows. In section 2 we review the model [1] and study the 
evolution of the homogeneous component of the infiaton. In section 3 we provide general 
formulae that allow to compute the spectra of scalar perturbations in this system. In section 
4 we give explicit expressions for the scalar power spectrum and bispectrum in terms of the 
main parameters of the theory. In section 5 we compute the spectrum of tensor modes and 
in section 6 we discuss our results. 

2 Natural steep inflation 

In order to make the paper self-contained, in this section we review the model and the results 
of \T]. We consider natural infiation with a pseudoscalar infiaton $ coupled to N' U{1) gauge 
fields. The Lagrangian density of the system is 

1 1 re ~ 

^ = -^{d^r - vi<^) - -F^,^Fr - ^$F;.Fr , (3) 

where V{^) = A^ [1 + cos($//)] with / <, Mp. The index a is understood to be summed 
upon, and it labels the ^(1) groups, ranging from 1 to A/". The parameter a is the dimen- 
sionless coupling constant between axion and gauge field. Introducing the vector potential 



A"(r, x) with a^ B" = V x A'*, a^ E" = —A'", we obtain the equations of motion 

d^ f 



( ^-V'-«-Vxj A'^ = 0, V-A" = 0, (4) 

where H = a'{T)/o?{T) and where the prime denotes differentiation with respect to the 
conformal time r. 

2.1 Generation of the gauge field 

The rolhng inflaton induces the generation of quanta of the gauge field. In order to study 
this process we promote the classical field A°(r, x) to an operator A"(r, x). We decompose 
A'* into annihilation and creation operators 

A^{t, x) = 5^ y — ^ [e,(k) AUr, k) a,(k) e^"^- + h.c] , (5) 



A==f 



where the helicity vectors e-i-(k) are defined in such a way that k ■ e± = 0, h x e± = ^i ke±. 
Then, the mode functions A^ must satisfy the equations A'^" + (fc^ =p a A; $'//) ^± — 0- 

Since we are looking for inflating solutions, we assume a{T) = —1/{Ht), and d^/dt = 
$0 = constant. Hence, the equation for A'^ reads 



d'Ai{T,k) , r 2 , .,e 



dr^ 
where we have defined 



+ 



k^±2k 



Al{T,k) = 0, (6) 



^ = "27^' (') 

and we will be interested in the case ^ ^ C^(l)- Depending on the sign of ^, one of the 
two solutions A" or A°i of eq. dol) will develop an instability. We assume that a > and 
V'{<^) < 0, so that ^ > 0. 

The solution that reduces to positive frequency for kr ^ —00 is 

Alir, A;) = -L= [zFo (i^, -kr) + Go (i^, -kr)] , (8) 

V2 A; 

where Fq and Go are the regular and irregular Coulomb wave functions respectively. At late 
times \kT\ ^ 2^ the positive helicity mode A^ behaves as 

]__ f ^ \ 7r^-2y/2 (k/aH 



^a / ^^ ~ " n^-2y/2^k/aH rQ\ 

and is thus amplified by a factor e'^^. On the other hand, the mode A°i is not amplified by 
the rolling inflaton and we will ignore it from now on. 



2.2 The slow roll solution 

The amplified gauge modes backreact on the rolhng infiaton, and slow down its rolling along 
V{^). Here we describe the portion of parameter space where such effect is so large to 
lead to slow roll inflation on a steep potential. In order to estimate the backreaction of the 
gauge field on the infiaton, we use the decomposition of A"(r, x) described above. Then, 
the quantity (E*^ ■ B*^) is very well approximated, for ^ ^ 3 (that, as we will see, is the regime 
we are interested in), by 

(E"-B'^) ^-2.4x lO-^AT ( — J e2"«, (10) 

that we plug into the first of eqs. (|4]), to obtain the equation of motion of $ (now written in 
the physical time t) 

^'* +3i7^ + n$) = -2.4 X 10-^ AT^ C^Xe^-^. (11) 



dt^ dt ' ' / V ^ 

Since we are interested in finding inflationary solutions where slow roll is supported by the 
dissipation into electromagnetic modes, we assume that both $ and 3if $ are negligible with 
respect to V'{^). In this case, the expression of ^ is controlled by the transcendental equation 



^ = i^'^ 



1.3x10^ M| /|r($)| 



(12) 



where we have assumed 3Mj,H^ = ^^^ + V{<^) + ^ (E^ + B^) ^ \/($). 

As we will see below, ^f = O(IO^), a = 0(10^) and ^ = 0(10), so that, using V ~ V/f, 
e-flogf^. 

As shown in [T], the slow roll conditions |<i>^ + ^ (E^ + B^) < V{^), \H\ < H"^, |$ + 
3H^\<t:\V'{^)\ can all be satisfied if 

When this condition is violated, the energy in the electromagnetic field cannot be neglected 
with respect to the energy in the infiaton. This happens when we are close to the bottom 
of the potential. Indeed, by approximating V{^) oc $^ near its minimum, we see that the 



condition (13) is violated when $ ~ i f /o.. This point marks the beginning of reheating, so 
that the reheating temperature is estimated to be Trh — A y/^Ja. 

One more condition on the parameter space of the system comes from the requirement 
to have enough inflation. The number of efoldings is given by 

*-f Hd^ a f^f d^ a $j - <l>i 



Since the range of variation of $ is bounded by |$j — $j| <, tt /, the above equation implies 
that a/^ > 2NJtt. 




Figure 1: Evolution of background quantities for A = 10"^ Mp, / = 0.1 Mp, a = 300 and 
A^ = 10^. Left panel: evolution of ,^(t). Right panel: the solid line correspond to the slow-roll 
parameter e = —H/H'^, the dashed line gives the ratio of the energy in gauge modes over 
the energy in the inflaton. In the inset we plot the number of efoldings as a function of time. 
The time t is in units of Mp/A^. 

How large can ^ be? Let us say that we require the reheating temperature to be larger 
than (9(10^) GeV, so to leave room for some mechanism of baryogenesis related to the 
electroweak phase transition. Then since ^/a ~ 10~^, we can take A as low as a few TeVs, 
corresponding to ^ ~ 20. In this case, since inflation occurs at a low energy scale, only ~ 30 
efoldings of inflation are needed to solve the problems of standard Big Bang cosmology, and 
we will need a ^ 400. In the opposite regime of very high energy inflation, ^ can be as small 
as 4 for A ~ 10^^ GeV - in which case A^e — 60 so that we need a ^ 150. As we will see in 
section 5, the requirement that the system does not lead to overproduction of tensors will 
impose stronger constraints on a. 



A numerical study of eq. (11) supports the validity of the analytical estimates above. 
We show in figure 1 the evolution of the scale factor, of ^ and of the slow-roll parameter 
e = -H/H'^ for the choice of parameters a = 300, A = 10"^ Mp and / = 0.1 Mp. The 
left panel of figure 1 shows that the parameter ^, rather than being constant, increases with 
time. The time dependence of ^ is however rather mild - it starts at .^ ~ 4 and ends at 
^ c^ Q after about 60 efoldings of inflation - and is in approximate agreement with the rough 
estimate ^ — - log — r^ that would give ^ ~ 4.4 for this choice of parameters. 



3 Perturbations 



Perturbations are usually generated as the quantum fluctuations of the inflaton are amplified 
by the evolving background. This is different from our scenario, where inhomogeneities in $ 
are sourced by those in the electromagnetic field, analogously to the situations encountered 
in warm [TOl [TT] or in trapped inflation [131 ES] . 

The curvature perturbation ^ on a uniform energy density hypersurface is related to the 
perturbation of the number of efoldings hj ( = 6^ = N{x) — N, where N is the number of 
efoldings in the homogeneous background. If we write the perturbed value of the inflaton 



field as $ = $o(t) + <P{t,x), then ( = H (p/^Q. The perturbation obeys the equation 

(P" + 2aH^'+ (-V' + a' V"{<^o)) = -y '^[E'^ " B"] , (15) 

where the fluctuation 6[Ei"' ■ B''] receives two contributions, the flrst from the intrinsic inho- 
mogeneities in E° ■ B" (that would be present for = 0), and the second from the fact that 
(E° • B") depends on <i>o 

,5[E'^ • B^ = [E'^ ■ B'^ - (E'^ ■ B'^)] ,^0 + ^(E°-B°) ^ ^ ^^^^ 

The above equation should also include a term proportional to through the dependence of 
ga . ga Qj^ jj^ Y)\ii this contribution turns out to be subdominant. 



We denote the term in brackets in eq. ( 16 ) by (5e.b, and the second term gives vr a V'(j)/{f H) 
upon using the background equation a (E*^ • B'') // = V . The Fourier transform of the per- 
turbation will then obey the equation 

For steep potential inflation we have it aV'/{f H^) ~ aMp/f^ » 1. Moreover, we are 
interested in the super- horizon wavelengths k « a VV". Hence, the solution of this equation 
can be written as 



(k, r) = -- / dn G(r, n) a\n) Jin, k) 

J J —oo 



the operator 



where J{t, k) = / '^-3/2 e *'''^(^e-b(t, x) , and the retarded Green's function G{t, t') for 

(27r) 






reads 



r 



z/ 1 — z/_ 



T \ '^+ / T 



r / \T 



GiT,T') = - - - - Qir-r'), (20) 



_ naV / A V'H^P 

^^ " 2 /if 2 (^1 ± Y 1 - ^ a2\/'2 

The second term under the square root in i>± scales as {f /aMpY and is much smaller than 
unity. As consequence 



aV M 



~a-^»l, 

«1, (21) 



fH^ P 

V"f 1 



Tt aV a 
where in the steps marked by the symbol ~ we have used V" ~ V / f ~ V/ p 
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The two point function of the metric perturbations has been computed, for our system, 
in pp. Here we will set up the calculations for the n-point function. See also [IHl [21] for a 
detailed analysis of this system in the regime where the backreaction of the gauge field on 
the inflaton is negligible. 

The n-point function of is given by 



(0(ki,r)0(k2,r) ... 0(k„, r)) = ( -- j dn dT2 . . . dr^ a'^ {ti) a"^ {T2) 
X Gir, n) Gir, r^) ... G{r, r„) ( J(ki, n) J(k2, r^) . 
Generalizing the results of yjj we obtain 

(0(ki,r)0(k2,r)...(/)(k„,r))=A/-aH (-J^ 5^^\k^ + k2 + ... + K) 



. c?{t^ 



(22) 



J(k„, r„)) 



X 



(i^q 



[e+(q) •e;(q-ki)I(r, |q|,|q-ki|)] 



(2vr)3"/2 

X [e+(q - ki) ■ e:;_(q - ki - k2)X(r, |q - ki|, |q - ki - k2|)] 
X---X [e+(q + k„)-e;(q)X(r,|q + k„|,|q|)] , (23) 



where 



X(r, /ci, /C2) 



dr 



'^^^' ^ {^1 A+(r', fci) AV(r', k2) + (fci O A;2)} 



a^fr' 



(24) 



and where o{n) = (2{n — l))!!/2" is a combinatorial factor. The correlation function (23) 
can be represented graphically as in Figure |2} 

Next, we compute X(r, fci, ^2)- Since we are interested in the spectrum at k << aH, 



we neglect the term {t /t'Y+ in the expression (21) of the Green's function. Then using the 
expression ^ for y4+, we obtain 



X(r, fci, k2) = f 



where 



se 



a(T)H 



{kl^^ + kT 



i72e2< k^l^ k, 



1/4 j,l/4 
2 



J 2"-+e (,1/2 + ^1/2 



{kr-(^.^/v.2^/f^-} 



,(25) 



/W 



yV^ 



(26) 



The quantity in curly brackets has a weak dependence on r, therefore we shall neglect its 
contribution to X. In addition, the argument of f[x\ can be set to zero as we are interested 
in correlations as r — )■ 0. In this case we have f[x\ = 7\. We will denote by X{ki, ^2) the 



function (25) with z/_ = and / = 7!, and we will use this expression in eq. (23) from now 
on. 



k + q — ki 

9 I q — k ki /^'^"~^ ^2 = —ki 



T 



(a) 




Q 
(b) 



t 9- ^^^9 + ^3 ^ 




(c) 



Figure 2: Graphical representation of the correlation (23). The vertex (a) denotes the 
quantity e+(q) • e;^(q — k)X(r, |q|, |q — k|) as a function of the momenta q and q — k. 
Diagrams (b), and (c) are the two-point and three-point functions. Dashed lines are the 
scalar perturbations 0(k). 



3.1 The two-point function 

The two-point function is 



(0(k,)0(k2)) = Ar^<5(3)(ki + k2)||-^[e+(q).e;(q-ki)X(|q|,|q-ki|)] 

X [e+(q-ki)-e;(q)X(|q-ki|,|q|)] . (27) 

In order to compute the scalar products of polarization vectors appearing in this equation, 
we must know the general expression of e+(q). Writing the components of q as 

q = g (sin 9q cos (f)q, sin 9q sin (f)q, cos 9q) , 

the polarization vector e-i-(q) reads 

. . 1 



Hence, we find 



V2 



(cos 9q COS 4>gT i sin 0„, cos 9g sin (t)q±i cos (/)„, — sin 9q 



(28) 
(29) 



|e+(q)-e;(q-k)| 



4 



1 + 



q ■ (q - k) 



n 2 



|q||q 



Computing numerically the angular integrals in eq. (27) we find 



(0(ki)0(k2)) = 2.1x 10-6 A/" 
that matches the result obtained in [Ij. 



a 



H^ e^"? 5(3)(ki + k2) 



p K e 



ki 



(30) 



(31) 



3.2 The three-point function 

The three-point function is given by 

X [e+(q-ki).e;(q + k3)X(|q-ki|,|q + k3|)] [e+(q + kg) ■ e;(q)X(|q + kg], |q|)] . (32) 



{<f>{kMk,)<f>{k,)) = -^ 5(=^)(ki + k2 + k3) I -^ [e+(q) • e;(q - k^ X(|q|, |q - k,\)] 



The correlator depends on the size and shape of the triangle formed by the vectors ki, 
k2, and kg. We denote 

|ki| = A;, \k2\=X2k, |k3|=X3fc, (33) 

with 

ki = (fc, 0, 0), k2 = -ki-k3, (34) 

k3 = -— (^1 - X2 + Xg, a/-(1 -X2 + Xs){l + X2- Xs){l - X2 - X3)(l +0:2 + X3), OJ . 

Using the equilateral configuration X2 = ^3 = 1, we find 

mo ^(k.) ^(k3)) ^ -lO-A-^ f g ^'""-' ;■-' + ■--' . (35) 

4 Computation of the power spectrum and f^i 

Given the equations above it is straightforward to derive the expression for the power spec- 
trum and the value of f^i for equilateral configurations in our system. The metric per- 
turbation C(k) is given by ({k) = H (j){k)/^Q. The scalar power spectrum can be read off 
from 

(C(ki) C(ki)) = V( p , (36) 



so that, using the background equation (12), the amplitude of the power spectrum turns out 

to be 

^ 5 X 10-2 

n-^i;^. (37) 

Comparing this result to the observed value P^ = 2.5 x 10^^ we find that JV ^'^ = 2 x 10^. 
Given that ^ cannot be larger than 20 or so, we deduce that we need Af ^ (9(10^) to explain 
the smallness of the cosmological perturbations in this scenario. 

Next, we can calculate f^l^ . This quantity is related to the three point function of the 
curvature perturbation by 

(C(kOC(k2)C(k3)) = ^ {2nf' f^,V'^S'-'\k, + k2 + k,)^-^, (38) 

that, using the formulae above in the equilateral limit ki = k2 = k^ = k, yields 

(C(k,)C(k.)C(k3)) = -0^3^^^!^('5i±^^^±M, (39) 
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and finally 



/. 



equil 

NL 



-1.3^ 



(40) 



equil I 



Since, depending on the energy scale of inflation, 4 < ,^ < 20, we have 5 < |/^" | < 30. In 
particular, if ^ — 20, as in the case of low scale inflation, then such a level of nongaussianities 
could be detectable by Planck. 

The scaling of the n— point function can be justified by the following argument (see 
also [21] for an analogous discussion of the two-point function in the case A^ = 1). Let us 
consider the main equation determining 0, eq. (17), in coordinate space and in cosmological 



time t rather than in conformal time r. The most relevant terms are the friction term, on the 
left hand side, and the source (5e.b on the right hand side. Considering that perturbations 
evolve for about a Hubble time before freezing out of the horizon, we can approximate 
~ if 0, so that we roughly obtain ^ (5e.b/^'- When we compute the correlators, since 
the various components of ^e-b sum incoherently, we can estimate (0") ^ (<^eb")/^'" — 
A/" (^e-b)"/^'"- Next we estimate ((5e.b) — (E ■ B) ~ / V /Ma, where we have used the 
background equation of motion given by the first of the eqs. (|4]). Finally, remembering that 
C = H (j)/^o, we obtain {("•) ~ A^ (^A/")""", that matches the dependence on ^ and M found 
above in the case of the two- and three-point function of (. It is worth noting that, in the 
language of |21] , the perturbations in this model are provided by a "feeder" mechanism that 
nonetheless leads to a hierarchical structure for the correlators. 



5 Gravitational waves 



The generation of gravitational radiation in this system has been studied in [19] in the regime 
where the backreaction of the gauge field on the infiaton was negligible. There it was found 
that the correlation functions of the left- and right-handed tensors take two contributions: on 
the top of the usual, parity symmetric part originating from the amplification of the quantum 
fluctuations of the graviton, there is a second, parity violating contribution generated by the 
presence of a classical gas of photons. The analysis of [19] applies also to the present situation, 
and gives 



V 



t,L 



H^ 



7r2M2 



V 



t,R 



H^ 



7r2Mp2 



1 + 8.6 X lO^W 



1 + 1.8 X 10-^A^ 



Ml e' 
Ml e 



(41) 



The tensor to scalar ratio r = {V^'^ + 'p^^R'^ j'p^ reads 

_ ^2 2 + 8.6xlO-Wi^fi 



7r2M|, 



Vc 



(42) 



By trading e^'"^ for the other quantities in the system using eq. (12), and remembering that 
pC ^ 0.05/(Ar^2), we obtain 



1 21/ 

^37r2M^ 



2.7 X 10' 



e fv'f 



a2 V V 



(43) 
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The first term on the right hand side of the above equation corresponds to the standard 
spectrum of gravitational waves originating from the fact that the inflating Universe is a quasi 
de Sitter space. The second term corresponds to the contribution to the tensor spectrum 
induced by the presence of a gas of photons inside the horizon. Since l^' / ~ \^, the current 
observational limit r <, 0.2 implies that V^l'^ ^ 2 x 10^^ GeV and a >, 35,^^. 

In the specific cases discussed at the end of section 2.2, this implies that for high scale 
inflation A ':^ 10^^ GeV, .^ ^ 4, we need a ^ 500 (larger than the requirement on a originating 
from the need of realizing a sufficient number of efoldings of inflation). In the case of low 
scale inflation A ~ 10^ GeV, one has ^ ~ 20 so that we need a ^ 10^. Again, this constraint 
on a is more stringent than that coming from the requirement of sufficient inflation. 



Let us note that the second term on the right hand side of eq. (43) - i.e., the dominant 
contribution to r when the gas of photons provides the main source of gravitons - can 
take a wide range of values depending on a and ^. This is different from the case with 
weak backreaction discussed in fiSl [T9] where r, for large ^, depends only on the slow- 
roll parameter. This discrepancy originates from the last term on the right hand side of 



eq. (16), that is negligible in the regime of weak backreaction, but becomes important in 
the situation studied in the present paper, and that affects the amplitude of the spectrum 
of scalar perturbations while leaving the tensor modes unchanged. 



Remarkably, eq. (43) shows that primordial tensor modes might be detected in the up- 
coming surveys (that are expected to increase the sensitivity to r by a factor of 10 or so) 
even if the scale of inflation is very low, provided the ratio ^^/a is large enough. In this case 
the tensors would be essentially fully chiral - a signature that would be observable in future 
surveys once primordial tensors are detected 



6 Conclusions and discussion 

The model proposed in [1] is able to produce inflation on a steep axion potential provided 
the coupling a is of order of a few hundred or larger. While in the case of a single family of 
f/(l) gauge fields the spectrum of scalar perturbations produced by this model is too large to 
account for observations, the presence of ~ 10^ ^(1) gauge fields can reduce the amplitude 
of the perturbations to the observed value. 

One notable property of this model is that, while it depends strongly on the parameters 
a and TV, it has only a weak dependence on the axion constant / and, more importantly, 
on the scale of inflation A. One might wonder whether the QCD axion could be a good 
infiaton candidate. We see two major obstructions against this idea. First, as we have seen, 
the reheating temperature Trh is of the order of A/10, that for a QCD axion would give 
^RH ~ 10 MeV. Whfle this is still large enough to allow for a successful Big Bang Nucle- 
osynthesis [22] , it is hard to understand how to generate the observed baryon asymmetry in 
the short time between reheating at ~ 10 MeV and Nucleosynthesis at ~ 1 MeV. Of course, 
one might imagine that baryogenesis happens during the final stages of infiation, with the 
gas of chiral photons providing a strongly parity- violating environment, one of the necessary 
conditions for baryogenesis. The strongest obstruction, however, comes from the fact that 
we need 10^ copies of the U{1) gauge field. By the time Nucleosynthesis occurs, all of these 
fields (with the exception of the one associated to the Standard Model photon) must have 
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decayed into ordinary matter. While this is definitely possible (if the gauge fields have a 
small mass, they can decay into matter) the short time window between T ~ 10 MeV and 
T ~ 1 MeV is probably too narrow to allow for this process to happen completely. 

In the present work we have studied the nongaussianities and the tensor modes generated 
in this scenario. The phenomenology of the scenario turns out to be rather rich: the (equi- 
lateral) bispectrum is consistent with the current limits, but its amplitude is sufficiently 
large to be detectable by Planck. In fact, the scenario described in this paper would be 
under severe pressure if future data were to constrain |/^^' | to be smaller than 5 or so. 
The non-observation of tensors imposes constraints on the parameter a that are about one 
order of magnitude stronger than those found in [1] . The detectability of tensors is strongly 
parameter-dependent. Remarkably, however, tensors might be large enough to be detectable 
(say with r >, 0.01 or so) even if infiation occurred at very low energies (similarly to other 
situations considered in |22]). As pointed out in j20j, this originates from the fact that, due 
to the conservation of angular momentum, vectors are (unlike scalars) an efficient source of 
tensors. More importantly, tensors produced in this scenario would be almost fully chiral 
- a property that might be detected in future surveys [26] and that would provide a very 
specific signature of this model, especially if chiral tensors were to be detected along with 
some degree of nongaussianity in the scalar spectrum. 
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